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ABSTRACT 

To understand better the early galactic production of Li, Be, and B by cos- 
mic ray spallation and fusion reactions, the dependence of these production 
rates on cosmic ray models and model parameters is examined. The sen- 
sitivity of elemental and isotopic production to the cosmic ray pathlength 
magnitude and energy dependence, source spectrum, spallation kinematics, 
and cross section uncertainties is studied. Changes in these model features, 
particularly those features related to confinement, are shown to alter the Be- 
and B-versus-Fe slopes from a naive quadratic relation. The implications of 
our results for the diffuse 7-ray background are examined, and the role of 
chemical evolution and its relation to our results is noted. It is also noted 
that the unmeasured high energy behavior of a + a fusion can lead to effects 
as large as a factor of 2 in the resultant yields. Future data should enable 
Population II Li, Be, and B abundances to constrain cosmic ray models for 
the early Galaxy. 



Subject Headings : Cosmic rays — Galaxy:evolution — nuclear reactions, nu- 
cleosynthesis, abundances 



1 Introduction 



The Population I abundances of 6 Li, Be, and B (LiBeB) have been thought 
for some time to have their origin in spallation and fusion processes between 
cosmic ray and interstellar medium (ISM) nuclei (see, e.g., Reeves, Folwer, & 
Hoyle 1970; Meneguzzi, Audouze, & Reeves 1971; Walker, Mathews, & Viola 
1985). In the past few years these elements have been sought in extreme 
Population II dwarfs, stars well known to exhibit the "Spite plateau" in 
lithium (Spite & Spite 1982), which is understood to indicate the primordial 
7 Li abundance (Walker, Steigman, Schramm, Olive, & Kang 1991). Recently 
these same Pop II stars have also been shown to contain beryllium, (Rebolo 
et al. 1988a; Ryan et al. 1990, 1992; Gilmore et al. 1992a, 1992b; Boesgaard 
& King 1993) boron, (Duncan, Lambert, & Lemke 1992) and most recently 
6 Li (Smith, Lambert, & Nissen 1992). These abundances provide important 
clues about the early galaxy, and cosmic rays have been considered the most 
likely production mechanism for LiBeB in these stars as well (Steigman & 
Walker 1992 (SW); Prantzos, Casse, & Vangioni-Flam 1992 (PCV); Walker 
et al. 1993 (WSSOF); Steigman et al. 1993 (SFOSW)). 

If indeed the extreme population II LiBeB arise from cosmic ray inter- 
actions, the study of their isotopic abundances opens an important window 
on astrophysics. In principle, we may be able to gain insight on early cos- 
mic rays, as well as early star formation rates and chemical evolution (e.g. 
PCV; Fields, Schramm, & Truran 1993; Silk & Schramm 1992). In addition, 
associated with these early cosmic ray events is an appreciable gamma-ray 
flux which would contribute (perhaps significantly) to the present diffuse 
gamma-ray background (Silk & Schramm 1992; Prantzos & Casse 1993). 
The early cosmic ray scenario is also of importance to cosmology. One may 
use the results of a galactic cosmic-ray (GCR) spallation and fusion model to 
help infer primordial 7 Li from the observed Li abundance (Olive & Schramm 
1993). One may also use the GCR model's cosmic ray spallation information 
on the Li isotopes, compared with their observed abundances, to deduce the 
amount of possible stellar depletion (SFOSW). 

Our previous work (WS, WSSOF, SFOSW) was an attempt at a relatively 
model-independent approach. Without assuming a specific model (of cosmic 
ray or galactic chemical evolution) we concentrated on testing the consistency 
of standard GCR models with the recent LiBeB observations and with pri- 
mordial nucleosynthesis where possible by examining elemental and isotopic 
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abundance ratios. In another approach, Gilmore et al. (1992b), Feltzing & 
Gustaffson (1994), Prantzos (1993), Prantzos & Casse (1993), Pagel (1993), 
and particularly PCV have chosen to adopt a particular detailed model of 
early galactic chemical and cosmic ray evolution to examine its predictions 
and compare with observations. 

In this paper we will examine in detail the dependence of the LiBeB 
abundances produced by GCR nucleosynthesis on the uncertainties of, or 
allowed variations in, the cosmic-ray model. As such, we will discuss the 
uncertainties in cosmic ray models. After a brief review of the data, we 
present some model options in section |3|. Our assumptions regarding the 
evolution of the abundances of a, C, N, and O are discussed in section |j. 
We address the issue of the Be and B slopes versus [Fe/H] in section^. The 
cosmic-ray spectrum and the confinement of cosmic rays will be discussed in 
section We also explore implications of these models on 7-ray production 
in section |7|. We draw conclusions in section 

2 LiBeB Abundance Data 

To make this work as self-contained as possible, we show in table I LiBeB 
isotopic abundances observed in Pop II halo dwarf stars. We list those stars 
in which at least two light elements have been observed as we will for the 
most part be primarily interested in elemental or, even better, isotopic ra- 
tios. The notation we employ below and use throughout the paper is that 
[X/H] represents the log abundance relative to the solar abundance, namely 
log(X/H) -log(X/H) Q and [X] = 12 + log(X/H). In the table, the iron abun- 
dance represents an unweighted "world" average. For the other abundances, 
a weighted average is given. The 6 Li abundance was taken from Smith et 
al. (1993). The 7 Li abundances were taken from Spite & Spite (1982,1986); 
Spite et al. (1984); Hobbs and Duncan (1987); Rebolo et al. (1988b); Hobbs 
& Thorburn (1991); and Pilachowski et al. (1993). The 9 Be abundances were 
taken from Rebolo et al. (1988a); Ryan et al. (1992); Gilmore et al. (1992a,b); 
Molaro, Castelli, & Pasquini (1993); and Boesgaard & King (1993). Finally, 
the boron abundances were taken from Duncan et al. (1992). 

The ratios of 6 Li to 7 Li, Li to Be and of B to Be are the observed ratios. 
Because, for Pop II, the dominant contribution to the 7 Li abundance comes 
from primordial nucleosynthesis rather than GCR nucleosynthesis, a certain 
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degree of caution is necessary when comparing the first two of these ratios 
to the predictions we discuss below. On the other hand because there is no 
appreciable big bang source for either Be or B (Thomas et al. 1993), the ratio 
of these two may be compared (unless there is an additional primary source 
for n B as discussed in Dearborn et al. (1988), Woosley et al. (1990), and 
Olive et al. (1993)). 

Table 1. Observed Pop II abundances of LiBeB isotopes 



Star 


[Fc/H] 


Li 


[Be] 


[B] 


b Li/ Y Li Li/Be* 


B/Bc* 


HD16031 


-1.9 


2.03 ±0.2 


-0.37 ±0.25 




251 ± 185 




HD19445 


-2.1 


2.07 ±0.07 


-0.14±0.1 


0.4 ±0.2 


162 ± 46 


3.5 ± 1.8 


HD84937 


-2.2 


2.11 ±0.07 


-0.85 ±0.19 




0.05 ±0.02 912 ±425 




HD94028 


-1.6 


2.10 ±0.08 


0.44 ±0.1 




46 ± 13 




HD132475 


-1.6 


2.05 ±0.09 


0.60 ±0.3 




28 ±20 




HD134169 


-1.2 


2.20 ±0.08 


0.71 ±0.11 




31 ± 10 




HD140283 


-2.6 


2.08 ±0.06 


-0.87 ±0.11 


-0.1 ±0.2 


891 ± 262 


6±3 


HD160617 


-1.9 


2.22 ±0.12 


-0.47 ±0.18 




490 ± 244 




HD189558 


-1.3 


2.01 ±0.12 


0.89 ±0.33 




13 ±11 




HD194598 


-1.4 


2.00 ±0.2 


0.37 ±0.12 




43 ±23 




HD201891 


-1.3 


1.97 ±0.07 


0.65 ±0.1 


1.7 ±0.4 


21 ±6 


11 ± 11 


BD23°3912 


-1.5 


2.37 ±0.08 


0.30 ±0.4 




117 ± 110 





*Ratios are extremely uncertain due to inadequate treatment of systematics. 
Formal errors on ratios are underestimates. 



Because of the importance of the B/Be ratio, we note that the values given 
in the table for Be and the ratio B/Be, are averages over observations by 
several groups. These observations themselves show some spread which may 
be significant. For the star HD19445, beryllium upper limits were obtained by 
Rebolo et al. (1988a), giving [Be] < 0.3, and B/Be > 1.3, Ryan et al. (1990) 
found the upper limit [Be] < —0.3 and hence B/Be > 5 which should be 
compared with the value of 3.5 in the table which represents the only positive 
identification of Be in this star by Boesgaard and King (1993) and is slightly 
discrepant with the upper limit of Ryan et al. For HD 140283, we have 
measurements of Be by three groups: [Be] = —1.25 ± 0.4 from Ryan et 
al. (1992) giving B/Be = 14 ± 14; [Be] = -0.97 ± 0.25 from Gilmore et 
al. (1992) giving B/Be = 7 ± 5; [Be] < -0.90 from Molaro, Castelli, & 
Pasquini (1993), giving B/Be > 6; and [Be] = —0.78 ± 0.14 from Boesgaard 
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and King (1993) giving B/Be = 5±3. Finally for HD201891, [Be] = 0.4±0.4 
from Rebolo et al. (1988a) giving B/Be = 20 ±26 and [Be] = 0.67 ±0.1 from 
Boesgaard and King (1993) giving B/Be = 11 ± 10. As one can see there 
appears to be a wide range in values (and uncertainties) in this ratio, and the 
values given my be completely dominated by systematics which are poorly 
accounted for in the stated error.Q In particular, ratios are most useful if 
constant surface temperature and surface gravity assumptions are used in 
the specific elemental determinations. Unfortunately, this is not yet the 
case. 

Indeed it is the uncertainties in these measurements and in the averages of 
these measurements which are themselves uncertain. Namely, the treatment 
of systematic errors is far from systematic. The conversion of line strengthes 
to abundances involves inferences on the surface temperature and surface 
gravity of the star. Most observational determinations have been made us- 
ing different sets of inputs. Though one can ascribe some uncertainty to 
chosen values of these inputs, it is not always clear to what extent these 
systematic errors have been incorporated into the quoted so-called statistical 
error, and different authors make divergent assumptions on the uncertainty 
of their assumed stellar parameters. Furthermore any average will surely un- 
derestimate the true error because of the mistreatment of systematic errors. 
Thus it is our feeling that B/Be and Li/Be ratios are extremely uncertain. 

3 Cosmic Ray Models for LiBeB Spallation 
Production 

Currently studied models for LiBeB production are based on the work of 
Reeves, Fowler, and Hoyle (1970) and the more detailed follow-up work of 

1 Since originally submitting this paper we have learned the results (Kiselman 1994) of 
a non- local thermodynamic equilibrium (NLTE) calculation of the boron abundance in 
HD 140283. This calculation suggests that the boron abundance calculated in the LTE 
approximation might underestimate the actual boron content in the star by by factors of 
3-4. Because this result relies on delicate atomic modelling, it should be regarded with 
caution; nevertheless, if this is the first indication of a high Pop II B/Be ratio, it can 
have intriguing implications for the sources of boron (see Fields, Olive, & Schramm 1994). 
In addition, the re-analysis and also points up the potential uncertainty in these difficult 
measurements. 
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Meneguzzi, Audouze, and Reeves (1971). They describe cosmic ray propaga- 
tion via the simple leaky box model. This assumes a spatially homogeneous 
distribution of sources, cosmic rays, and interstellar material. Within this 
model, the propagation equation is 

dt dT T e ff 

Here Na = A^(t, T) is the number density of cosmic ray isotopes A = LiBeB 
at time t with energy (per nucleon) between T and T + dT. Note that eq. (|l|) 
is the general propagation expression for A, a primary or secondary species, 
with different terms being important for each species type. The first term 
on the righthand side, J a, includes all sources for A: 

Ja{T) = Q A (T) +5> / dT'fa(T')-£ (T, T') . (2) 

IJ th 

Here Qa(T) is a possible galactic source of A, given as a number rate per 
volume and per unit energy; and spallation production of A appears in the 
sum which runs over projectiles i and targets j. Additionally, fa = NiVi is 
the cosmic ray flux spectrum of species i and, af- is the cross section for 
the process i + j — *■ A + ■ ■ •. The second term in the propagation equation 
(|TJ) accounts for energy losses to the ISM, with b A = —(dT/dt)A allows for 
ionization from the Galaxy. The third term of the propagation equation 
accounts for catastrophic cosmic ray losses, 

:r- = E + (3) 

T e ff i 'esc 

with a^ 1 encoding spallation losses of A in the ISM, and r esc (T) being the 
lifetime for cosmic rays against escape. 

The propagation equation ([!]) is solved for the case of a steady state, 
dN/dt = 0, in which cosmic ray production is in equilibrium with the losses. 
One thus obtains the spectrum N of these elements, propagated from their 
source J. At present, we will ignore losses due to inelastic nuclear collisions, 
valid when r~l 3> ni<j™ A el VA (but see section |6TT| ). We assume that the pri- 

and CNO, have some homogeneous galactic 
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source J = Q, and negligible spallation production or losses: a = 0. Writing 
the solution for these in terms of the cosmic ray flux = Nv, we have 



1 f°° 

MT) = 777a / dT' qi (T') exp (- [R t (T') - it> (T)] /A) (4) 
Wi(T) Jt 

where W{ = h/p IS MV, and g, = Qi/pisM, and i — p, a, CNO. Also, Ra(T a ) = 
Jq A dT' / (dT' /3X)ism is the ionization range which characterizes the average 
amount of material a particle with energy T A can travel before ionization 
losses will stop it (expressed in gem -2 , as X = pisAivt). The ionization 
ranges were taken from Northcliffe and Schilling (1970) for low LiBeB energies 
(< 12 MeV/nucleon) at which partial charge screening effects are important, 
and from Janni (1982) for higher energies, using the scaling law 

R A (Z;T) = A/Z 2 R P (T) . (5) 

Here T is the kinetic energy per nucleon and R p is the proton range. 

In figure [1] we plot the proton and a fluxes <pi calculated from eq. ((4|) . To 
show the effect of energy losses on the propagated flux, we also plot Ag^, the 
solution of eq. (0) for negligible energy losses {pi ^0). As we will discuss 
below in sections [5] and 6.1 , these losses are important at low energies and 



negligible at high ones. This is manifest in figure [I], which shows the scaling 
= Ag to be followed closely at high energies, while at low energies the 
ionization energy losses significantly reduce the propagated flux from this 
scaling. This behavior is qualitatively similar for the two source types we 
consider; propagation differences at low energies are discussed in section 164] . 

In solving the propagation equation for the secondary elements LiBeB, the 
source term J is assumed to have no primary source component, i.e. Q = 0. 
The exclusive production of these elements occurs via spallation between the 
primaries and the ISM: of^ ^ 0. One can write an expression like eq. (||) for 
LiBeB, giving their cosmic ray spectrum. What we wish to know, however, 
is not this steady-state spectrum but instead the amount of LiBeB thermal- 
ized and added to the ISM. To compute this from the LiBeB spectrum one 
assumes that all such nuclei below some threshold kinetic energy T t herm are 
thermalized and then one examines the LiBeB current ^^^(Tiherm) below 
this threshold. Below the lowest spallation threshold there is no source term 
in eq. ([!]) for LiBeB, and the ionization loss term is much larger than the 
escape term. Thus to a good approximation the propagation equation ([I]) 
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reads 



d_ 
df 



b A N A (T)=0 ,T<T t 



th 



mac 



(6) 



and so the LiBeB current ^^(^therm) is constant for T t h er m below spallation 
production thresholds. One may thus choose any T t h erm < T™ c at which to 
evaluate this current; we choose ours to be right at threshold. 

With this method of computing the production rate of LiBeB via eval- 
uation of the subthreshold LiBeB current, we can write the rate of LiBeB 
accumulation in the ISM as 



where A = 6,7 Li, 9 Be, 10 B, i/a = n A /n>H, the 4>i come from eq. (f|), and we 
have ignored the small time variation of nu- 

The factor Sa{Ta, t) accounts for the energy loss of nucleus A in the ISM 
and gives the probability of its capture and thermalization; it is a function 
of the lab energy T A of the daughter nucleus A, as well as the epoch t, and 
is given by 



These energy losses in the ISM serve to trap the spallation products and so 
to allow them to contribute to the LiBeB finally thermalized in the ISM. This 
process competes with the the spallation products' escape from the galaxy, 
quantified by the escape rate r~l of eq. ([[]). This quantity appears in eq. (§) 
through A = pisMVT esc the (possibly) energy- dependent average pathlengthQ, 
also in gem -2 . Note that as r esc sets the scale for the cosmic ray residence 
time before escape, A is the amount of matter traversed before escape. As 
any T therm < T^ c gives R(T therm ) < A, we may put S A ^ exp(-R A /A). 
For "forward" kinematics, with a light cosmic ray nucleus impinging on a 
stationary, heavy ISM nucleus, Ta is small for typical spallation energies, 
and so Ra <C A and S A — 1- For "inverse" kinematics, with a heavy cosmic 
ray nucleus on a light ISM particle, S A can differ significantly from unity and 
the LiBeB yields are reduced accordingly. 

2 By writing the argument of the exponential in Sa as we have i n eq . (JsJ) , we are 
assuming A to be constant in energy. For further discussion see section 3.1 . 




S A (T A , t) = exp [- {R A (T A ) - R A (T therm )} /A(t)] 



(8) 
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Following Walker et al. (1992), we put 

qi (T,t)=yf R (t)q p (T,t) (9) 

with % = p,a, CNO: i.e., we posit the constancy of the cosmic ray isotopic 
and elemental ratios at the source and over the entire energy spectrum. We 
then choose to make the more serious assumption that we may express the 
proton source strength in the separable form 

q p (T,t) = q p (T)f(t). (10) 

We will in fact consider various Population II source spectra q (see section 
|6.2j ), but we will not allow a fully general energy or time dependence. While 
this has the immediate advantage of simplifying the calculations, one may 
view this assumption as a postulation that the present mechanism of cosmic 
ray acceleration does not differ dramatically over time in its energetics, but 
only in its net cosmic ray output. 

In our analysis of cosmic ray model features, and in the accompanying 
figures, we will concern ourselves with the effect on the LiBeB production 
rates — or rather their ratios — as calculated by eq. (0). Given a set of cosmic 
ray and ISM abundances, and a confinement parameter A, these rates may 
be evaluated numerically. As in previous work, we used empirical values of 
partial cross sections tabulated in by Read and Viola (1984). To proceed 
further and integrate the rates to get the LiBeB yields would require a full 
model for cosmic ray and chemical evolution. We discuss this issue in the 
next section. 

As a cosmic ray model of propagation, the above amounts to the simple 
leaky box (reviewed by, e.g., Ceasarsky 1980, 1987). This is the simplest 
model used to describe cosmic rays today, and as such has obvious compu- 
tational advantages. While the simplicity of the leaky box makes it a useful 
tool, it is physically unrealistic and thus, in some cases, quite inaccurate. A 
proper treatment of cosmic ray propagation must explicitly include diffusion 
effects only sketched with the leaky box confinement parameter A. Also, a 
proper model must abandon the simple leaky box assumption of spatial ho- 
mogeneity of sources and interstellar material. Despite these shortcomings, 
however, we will follow previous authors and adopt this model in what fol- 
lows, both for the above pragmatic reasons and moreover because the epoch 
we consider is so poorly understood in its relevant details that adoption of a 
more detailed description is not warranted at the present time. 
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4 The Role of Chemical Evolution 



Although we can calculate the rates in eq. (|^), for a given set of parameters, 
using well-understood physics, to integrate these rates requires knowledge of 
the chemical and dynamical evolution of the cosmic ray and ISM species in 
the Galaxy's early history. Such knowledge is sketchy at present, but some 
reasonable assumptions may prove fruitful. 

If we assume that the H and He abundances do not change much from 
their values as set by the big bang, we have 

H{t) « H BB 
y He (t) « yl B « 0.08 . (11) 

If we further assume that the CNO nuclides evolve at the same rate, so that 
their ratios remain constant within the early epoch considered here, we put 

Vcit) _ y N (t) _ ypjt) 

„,OBS „,OBS „,OBS ■ y LZ > 

Dc Vn Vo 

Here we denote with OBS the observed CNO abundances for a given Pop 
II Be or B star. Equation (|P2| ) is assumed to hold for times t less than the 
galactic age r at the birth of the star. One should interpret equation [12] 
with some care, as approximation of constant C:N:0 ratios is not strictly 
true. We expect these elements to have different sources: O is made in type 
II supernovae but C and N are made primarily in intermediate mass stars. 
Thus any differences in the evolution of these sites will result in differences 
in the abundance ratios. Fortunately, large differences between CN and O 
production are only important at early times, and the observed Pop II CNO 
abundances are consistent with constant C:N:0 over the timescales important 
here. One should note, however, that in Pop II stars the O/CNFe ratio, 
though roughly constant, does exceed the solar ratio by a factor of ~ 3 . This 
comes about because of the predominance at this epoch of type II supernovae 
over type I's. The effect is important for our purposes, as spallation of O 
gives a lower B/Be ratio than spallation of C. 

We again simplify by extrapolating the present relation between abun- 
dances of cosmic ray and ISM species, namely we assume that the relative 
enhancement of the cosmic ray CNO abundances over interstellar CNO abun- 
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dances has remained constant in time. We put 

SH?) < i3 > 

" 2 V / \ & % / present 

where i =CNO and the righthand side denotes the present value. We will in 
fact adopt the stronger assumption, made in most works, that 

y CR (t) = y ISM {t) (14) 

but we note that this assertion is false for the present cosmic rays. In partic- 
ular, Asakimori et al. (1993), and references therein, report a H/O depletion 
in present cosmic rays of about a factor of 2 relative to the solar system 
ratio. The proton-to-helium ratio in the cosmic rays is energy dependent 
and varies from ~ 0.2 at the lowest measured energies to ~ 0.05 at the high- 
est energies which interest us here. Additionally, as shown in Buckeley et 
al. (1993), He/O in the cosmic rays is depleted by a factor of ~ 4.6 relative 
to its local galactic value. In this paper we have normalized to the protons, 
and thus adopting the observational results would amount to reducing ?/ H e 
by 2.3, and enhancing ?/cno by a factor of 2. We will instead use eq. fly}) 
to allow comparison with previous results. Bear in mind, however, that the 
cosmic ray abundance scalings are uncertain by at least a factor of 2. 

There are proposed scenarios in which early cosmic ray abundances do 
not scale with early ISM abundances. The two suggestions are constructed to 
allow for spallation to occur with CNO abundances at high, Pop I levels, thus 
providing a natural explanation for the Be/O constancy. The first of these 
posits that supernovae accelerate their own (metal-enriched) ejecta, which 
makes the cosmic ray CNO abundance high. In the case, LiBeB production 
is dominated by the spallation of these fast, heavy nuclei, rather than by 
spallation of fast protons and a particles. We may model this behavior by 
putting 

yCR{t) = ysupernova „ present > y OBS ^ 

The other suggestion is that spallation occurs predominantly in (temporarily) 
metal-enriched regions around supernova remnants. In this case the cosmic 
ray abundances do scale with those of the ISM, since the target abundaces 
are always at Pop I levels. In what follows, however, wherever not explicitly 
noted otherwise, eq. (TT3|) will be assumed to hold. 
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If we allow the confinement pathlength to vary with time, then as we will 
see below, we may not further simplify the rate equation (|^), or the integral 
that is its solution. If the confinement is taken to be constant in time — a 
dubious proposition, as we shall see — we may now separate the propagated 
flux: 

<j> i {T,t)=y? R f{t) V {T) , (16) 

where yf R is the cosmic ray abundance of i relative to protons, as in eq. @. 
This allows one to write the solution to the rate equation as a sum of terms, 
each of which has a factor that involves an integral over time and a factor 
involving an integral over energy: 

y A (r) = £ *V ( 17 ) 

ij 

with the "exposure time" given by 

A tJ (r)= I" dt yj (t)yf R (t)f(t) (18) 
Jo 

with r is the age of the Galaxy at the birth of the star in question, and the 
"reduced rate" is 

= / dT <fi(T) af^T) S A (T A ) . (19) 

th 

Note that these factors are different for forward and inverse kinematics, that 
is, ($cr^)/ 7^ ($<7^) r . As mentioned above, this stems from the longer range 
and greater chance of escape for the fast A-nuclei produced in inverse kine- 
matic collisions. 

With the assumptions made thus far in this section — most importantly, 
that of the time constancy of A, the number of independent exposure times 
may be reduced to three (SW, WSSOF): 

i = j = a 

Ap^ f forward kinematics (20) 
ApQ inverse kinematics 

Of these, we will take A QjCf as the fiducial quantity, as it measures the in- 
tegrated flux enhancement exclusively (i.e. with no explicit dependence on 
chemical evolution), and we define 

T ^ A a Jy 2 He = f T dt fit) (21) 

J 
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We then define the quantity (0/H), which is a measure of the average abun- 
dance of oxygen relative to hydrogen,: 



(0/H) ISM ee A>% M /T 



p,o 

£ dt yff M (t)/(t) 
Jo dt fit) 



(22) 



with a similar expression for (0/H)cr- Note that (O/H) is to the weighted 
average of O/H at time r. This average weighted by the net flux enhancement 
f{t), and is normalized. 

With this notation we have 

y Be ,B = TY,'y? R y] SM (23) 



X 



I ,,/SM I J 1 CR 

\ yo yo / 

m = r [yL<*<&> (24) 

ii \ Vo Vo J 

where the prime on the summation indicates that any i = j = a term has 
been deleted. Note that the ratio of LiBeB production rates, i.e. of eqs. (p4]) 
and (j25D , does not depend explicitly on A aa . However, the flux enhancement 
appears via the strong dependence on (O/H), which includes a factor of A aa 
in its definition. 

From eq. fiTTD it is clear that chemical evolution effects are encoded into 
the LiBeB abundances in terms of exposure time, in which these effects 
appear in convolution with cosmic ray flux evolution. While this means that 
the absolute abundances are dependent on chemical evolution, it also means 
that the relative abundances are much less so if the LiBeB isotopes compared 
originate in similar processes. This is the case with Be and B, which are both 
produced in pure spallation events, as opposed to Li, which is produced not 
only by spallation but also in a + a fusion events. This point was stressed 
by WSSOF, who noted that the ratio B/Be is the least dependent on our 
knowledge of the early galaxy and of early cosmic rays. WSSOF showed that 
any model for early cosmic ray LiBeB production yields a B/Be ratio lying 
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within the range 7.5 < B/Be < 17. One might compare this with the average 
of the observed values in the most extreme Population II dwarf, HD 140283, 
which has B/Be = 6 ± 3. This ratio is the strongest prediction of the cosmic 
ray model. However, while the theoretical prediction is firm, particularly in 
setting a lower bound to the B/Be ratio from cosmic rays, the experimental 
ratio is still difficult to obtain reliably. As discussed in section |2], the true 
errors in the measured value, when including systematic uncertainties, are 
certainly larger than the simple weighted mean we quote here. Comparisons 
with the observations are thus not all that compelling at the present time. 

Under the assumption that there is no additional (primary) source for 
n B, the B/Be ratio is weakly dependent on chemical and cosmic ray evolu- 
tion. However, it is well known that the predicted value of the 10 B to 11 B 
ratio in GCR nucleosynthesis is about 2.5, whereas the observed ratio in 
meteorites is just over 4. Dearborn et al. (1988) proposed a possible ther- 
monuclear production of n B in Type II supernovae, and Woosley et al. (1992) 
had suggested that neutrino processes during supernovae may yield a signif- 
icant abundance of n B. Both suggestions are tested in a GCR and chemical 
evolution model by Olive et al. (1993), where it was found that while the 
11 B/ 10 B ratio could be brought into agreement with the observation, the su- 
pernova sources added a significant primary source for n B which upsets the 
constancy of the B to Be ratio as a function of metallicity. Though the data 
seem to disfavor this primary source which predicts a B /Be ratio in excess of 
50 at [Fe/H] < —3, the scant amount of data at present can not conclusively 
eliminate this possibility. 

In contrast, the ratios of Li with Be and B are always strongly dependent 
upon chemical and cosmic ray evolution because of the a + a channel for 
production of Li. This is in fact the most important process for Pop II 
cosmic ray synthesis of Li, since in the metal-poor Pop II environment, O/H 
can be reduced from its solar value by two or more orders of magnitude, 
while He/H is undepleted. In the language of eq. (O/H) <C y 2 He) and 
so the only important contribution to Li synthesis comes from the a + a 
term. In addition, the dominant component of the observed 7 Li is primordial, 
further complicating comparison of Pop II LiBeB with GCR nucleosynthesis 
predictions. 
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5 The Be and B versus Fe Relations 



Many have suggested (e.g. Vangioni-Flam et al. 1990) that since Be and B 
arise from the interaction of stellar nucleosynthesis products with cosmic 
rays, they should vary quadratically with elements, e.g. O, which form their 
targets. (Recall that while O/CNFe is enchanced in Pop II compared to 
Pop I, the ratios are still roughly constant and so a relation between LiBeB 
and CNO also manifests itself as the same relation between LiBeB and the 
more commonly observed Fe.) This quadratic dependence is argued to follow 
from the fact that cosmic ray sources are presumed to be supernovae which 
are the producers of the oxygen targets. However, cosmic rays have a large 
confinement time and hence the instantaneous recycling approximation that 
yields this quadratic dependence is quite inappropriate. Even more signifi- 
cant is the fact that the data does not show this dependence and is instead 
best fit by a Be slope relative to O that is much closer to 1 than to 2. This 
discrepancy, has received some attention and has even been interpreted as a 
failure of the cosmic ray model. Again, this conclusion ignores the effects of 
cosmic ray confinement, as we now demonstrate. 

Let us now go through the logic in greater detail. The argument for a 
quadratic scaling of Be and B with O proceeds from these assumptions: (a) 
Type II supernovae are the dominant source for O, so one expects a relation 
to Nsn, the total integrated number of supernovae, of the form O oc N$n- 
(b) Supernovae are also likely to be the acceleration site for cosmic rays, so 
one might naively expect the cosmic ray flux to vary with the current rate 
at which these acceleration sites are created: oc dN SN /dt ) if we ignore 
confinement and assume instantaneous recycling of the cosmic rays. Put 
slightly differently, in the leaky box model we have adopted, we have assumed 
an equilibrium to hold between the rate at which cosmic rays are lost and 
the rate at which they are created. In this model one can sensibly argue that 
the cosmic ray flux is indeed proportional to this creation rate. Then, taking 
the essential pieces of the rate equation ([?[), we obtain the naive result 

dBe/dt oc 00 

oc NsNd NsN/dt 
=4 Be oc O 2 (25) 

which gives the advertised quadratic dependence, with an analogous equation 
for B evolution. 
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While assumption (a) seems solid, assumption (b) is not because it ac- 
counts for evolution of cosmic ray sources but not for changes in confinement. 
Indeed the evolution of A is the critical feature allowing one to avoid the naive 
secondary-to-primary ratio for the B- or Be-to-0 ratio of eq. For a given 
source density q of primary cosmic ray elements (at some instant t), the prop- 
agated flux (f) increases with A, as at larger A the cosmic rays traverse more 
matter before escape and so the steady state between production and escape 
is reached at a higher flux. In the limit of no energy losses, valid at high 
energy, the relation is = Aq. With energy losses, this relation is energy 
dependent, but for the total flux PCV find that it varies roughly as oc A 1 / 2 . 
Thus we should alter the cosmic ray flux scaling to be <f>(t) oc A 1 / 2 (t)NsN(t), 
and it follows that changes in the early confinement alter the cosmic ray flux, 



and its associated LiBeB production. Specifically, in the place of eq. £5|, we 
have the more general relation 

dBe/dt oc A n OdO/dt 
^dBe/dO oc A"0 (26) 

with n < 1 the index for the scaling law between cosmic ray flux and con- 
finement (with PCV finding n = 1/2). If A can vary with time — and indeed 
it is hard to imagine how to avoid such variance — then the simple quadratic 
relation between Be or B and O is clearly not the prediction of the cosmic 
ray model. The effect of the variance of A is further discussed in the next 
section. 



6 Cosmic Ray Model Dependences 

We wish to investigate the dependence of calculated elemental ratios Li/Be 
and B/Be on various model parameters, for different epochs, i.e. differ- 
ent [Fe/H]. Such a calculation, however, requires specific models for galactic 
chemical and cosmic ray evolution. Without constructing such models, what 
we can examine, at a given [Fe/H], are the ratios of production rates 

pl^^dU/dBe and = dB/dBe (27) 

dy B e/dt dy Be /dt 

(c.f. eq. (0)) and their dependence on escape length. 
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This comparison of production rates provides useful bounds on the ele- 
mental ratios any full model would predict. Specifically, at a given epoch 
[Fe/H] and confinement A, dLi/dBe gives the value Li/Be would take if the 
ISM abundances and cosmic ray confinement had been constant for the en- 
tire history of the galaxy up to that epoch. The ISM abundances of heavy 
elements CNO certainly are not constant but build up from zero, while the 
ISM abundances of H and He do remain roughly constant. The chemical evo- 
lution of the ISM affects Li and Be or B differently; Li is largely unaffected, 
being made mostly by a + a fusion, while Be and B are made by spallation of 
the CNO targets whose abundances will have always been rising to the levels 
observed at a particular [Fe/H]. The calculated dLi/dBe will thus provide 
a lower bound on the actual Li/Be a full chemical evolution model would 
predict for a given [Fe/H].0 The B/Be ratio, however, does not change at 
all with CNO evolution as long as the C:N:0 ratios remain constant, and 
thus is equal to dB/dBe. These ratios are roughly constant over the Pop II 
epoch we consider here, but they do of course change eventually as the O 
overabundance diminishes to it Pop I value. 

6.1 Confinement 

PCV have detailed a chemical evolution model that can fit observed Be and 
B data and does produce a slope versus iron that is smaller than the naive 
quadratic dependence. In their model, the pathlength is assumed to be larger 
than the present ~ lOgcm -2 , due to the larger confinement volume of the 
halo phase versus the disk phase, as well as the larger early gas fraction. 
They suggested that the confinement varies over the range 10 g cm -2 < A < 
1000 gem -2 . This variation arises specifically from scaling the escape length 
as A oc a g H; where the gas mass fraction a g decreases with time according 
to the adopted chemical evolution, while the disk scale height H has an 

3 Just as the CNO abundances will have changed, so too the cosmic ray flux strength 
and confinement will have evolved. Though an overall change in flux strength does not 
affect elemental ratios, evolution in confinement will change the spectrum, with high early 
confinement leading to a harder spectrum in the past (see next section). This decreases 
the Li/Be ratio, as the a + a production falls off at high energies. This goes in the opposite 
direction as the chemical evolution effect. However, we will in general explicitly examine 
the allowed confinement parameters, and in this context we may understand the dLi/dBe 
ratio to be bounded from below. 
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exponentially decreasing evolution imposed on it by hand, to simulate the 
dynamics of the halo collapse. 

While this line of reasoning is attractive, one must be cautiously aware 
of the uncertainties involved. For example, the early behavior of the main 
confinement mechanism, the galactic magnetic field, is not well understood. 
The favored means for field generation is through dynamo action tied to the 
galactic rotation, setting a timescale for generation of r rot ~ 10 8 yr. Note 
that in this case the galaxy does not as efficiently retain LiBeB long enough 
to thermalize them, particularly given the flatter (harder) cosmic ray spectra 
and/or the inverse kinematics some have suggested (see section |6.2| ). 



Furthermore, Malaney and Butler (1993) have pointed out that the escape 
length is bounded from above. Physically, this happens when the confinement 
becomes so large that nuclear destruction of fast LiBeB becomes important. 
In terms of eq. ([TJ), we have 

M + TT < 28 > 

where A~^ c = J2i n %^TA I Pism- Malaney and Butler note that A nuc is species- 
dependent (whereas A con f is not), and they estimate A nuc < 200 gem -2 for 
protons, thus setting an upper limit for A. Prantzos and Casse (1993) have 
included this effect in the PCV model and find it makes little change in their 
results. 

What is the effect of a changing confinement? A variation in confinement 
acts not only to alter the flux intensity but also to change the shape of 
the spectrum. This behavior is apparent in figures ||(a) and (b), which plot 
spectra for a fixes source strength q but with different values of confinement 
A. Notice that in both plots the total flux is higher for larger confinement, 
but by different amounts in different energy regimes. The lowest energy 
(< 100 MeV/nucleon) portion does not enjoy the same enhancement as the 
high energy (> 1 GeV/nucleon) region. One can see in eq. (f|) that the 
spectrum turns over roughly at the energy at which R(T) ~ A. A spectral 
peak at a larger T means increased power at higer energies. In other words, 
increasing confinement enhances the highest energy regime and so hardens 
the spectrum. 

Note the different low-energy flux behaviors for the two source types (fig. 
HI), the forms of which are given explicitly below (section |6.2| , eqs. (|32|) and 
Q33|)). The momentum source (broken curve in fig. |l|(a)), diverges as T — > 0, 
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whereas the energy source (broken curve in fig. []](&)), is finite for small T. 
These differences manifest themselves in the propagated fluxes. The flux 
(solid curve in fig. 0(a)), derived from the momentum source has a mild 
turnover at small energies, and the different A fluxes all converge. This 
behavior is contained in the flux equation (|J). At low energy (T <C m p ), 
the escape length R{T) ~ T/(dT/dX) <C A, is slowly varying, and so the 
exponential in eq. (|j) changes slowly as well. On the other hand, the source 
is a featureless and steeply falling power law: q{p) ~ p~ 2 ~ Thus one 
can, at sufficiently low energy, ignore the confinement parameter altogether 
and obtain 

^"^F)Sr dr ^ (29) 

which is independent of confinement. For a source power law in total energy, 
however, the source has a characteristic scale m p and does not strongly vary 
over this scale. Since R(T)/A does change appreciably in this range, eq. ([29]) 
is not appropriate, and one must use the full flux equation, with its the 
strong dependence on confinement. Thus the propagated fluxes (solid curve 
in fig. p](fr)), do not display the same convergence as those derived from a 
momentum-law source. 

The well-known (e.g. PCV) dependence of the propagated flux on con- 
finement has important consequences for elemental ratios. Although changes 
in the total cosmic ray intensity will not affect elemental production ratios, 
spectral changes can. One calculates LiBeB production rates by integrating 
the product of cross section with flux over kinetic energy. Thus the energy 
regimes in which this product is the largest will contribute most to the pro- 
duction rates. As shown if fig. [L], the flux peaks either at ~ 100 MeV/nucleon 
(for a momentum spectrum), or at ~ 1000 MeV/nucleon (for a total energy 
spectrum). This behavior is to be cmopared to that of spallation and fusion 
cross sections, which typically display low energy resonance peaks up to a 
few hundred MeV/nucleon, followed by a high energy asymptotic behavior 
characteristic of the interaction. Specifically, the cross sections approach a 
constant at high energies for spallation interactions with a many body fi- 
nal state, but suffer an exponential dropoff for fusion reactions which have 
few-body final states. (It should be cautioned here that the continued ex- 
ponential falloff of the a + a fusion reactions above a few hundred MeV 
is only an extrapolation since no data exists at higher energies; c.f. section 
|6.3| .) Thus a change in confinement has the effect on dB/dBe — the ratio of 
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pure spallation products — of changing the emphasis on the resonance region 
relative to the plateau. Both p and a cross sections for B production have 
lower thresholds than those for Be, and B cross sections also display higher 
peaks relative to their plateaus. Consequently, the dB/dBe ratio is lowest 
at cosmic ray energies in the plateau region. To produce a low B/Be ratio 
requires a spectrum that gives a larger weight to this energy regime, i.e. a 
harder spectrum. Just such a spectrum arises from a large confinement scale. 
The absolute effect on the dB/dBe is small, but important given that B/Be 
is constrained theoretically and observationally to lie in a small range. 

For the dLi/dBe, B production ratio, spectral flattening via large confine- 
ment leads to more pronounced changes. The effect again stems from the 
high energy cross section behavior: while the Be and B cross sections are 
asymptotically flat, the a + a — > Li fusion cross sections (which dominate 
the Pop II cosmic ray production of Li) are assumed to fall off exponentially 
at high energy. Thus harder spectra have a strong effect dLi/dBe as they 
emphasize an energy regime where the Li production is not just reduced rel- 
ative to Be, but in fact negligible. Thus dLi/dBe strongly decreases with 
increasing confinement. 

We have computed numerically the ratios of rates for the range 10 g cm~ 2 < 
A < 1000 gem -2 . The results appear in figures and|4], and are in agreement 
with our expectations. The strong dependence of these ratios on the differ- 
ent forms of source spectra will be further discussed below (section |6.2|) . In 
reading these figures and comparing to observational data, one must bear in 
mind the need to include evolutionary effects properly in order to go from the 
plotted ratios of rates to elemental ratios. Furthermore, as argued above, we 
see the falloff of both ofLi/effieand dB/dBe at large A. Note that the scales 
in these figures are very different. As A is increased from 10 to 1000 gem -2 , 
dLi/dBe falls off by a factor of 10 (fig. |3]), whereas dB/dBe drops off by 
only about 15% for the momentum source spectrum and even less so for the 
energy source spectrum (fig. f|). 

In figure [5] we plot the production ratio dLi/dBe against [Fe/H]. Such a 
plot suggests a comparison with observational data, but bearing in mind the 
caveat that the dLi/dBe value is only a lower limit to the actual value a full 
chemical evolution model would give. See the end section 6^2 for a discussion 
of such a comparison. Also note that a plot of dB/dBe as a function of [Fe/H] 
is uninteresting as long as one assumes the C:N:0:Fe ratios remain constant 
during the Pop II epoch, since in this case the dB/dBe ratio remains constant 
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(for a fixed confinement strength). 

In our considerations thus far we have followed most previous work on 
early LiBeB in putting A to be a constant in energy. This parameterization 
offers a simple way of trying to characterize the early galactic conditions 
which, as has been pointed out, are poorly understood. While offering sim- 
plicity, this approach does suffer the criticism that the choice of a constant 
escape length does not fit the observed present day cosmic ray pathlength 
distribution. As shown most recently in Garcia- Munoz et al. (1987), the 
escape length varies with energy, increasing slowly to a maximum around a 
GeV/nucleon, then decreasing as T~ 06 . They find that their empirical form 
for A(T) is well fit by 



and j3 = —0.6 for T > T\. A plot of this function appears in fig. |]. As 
is well-known, an energy-dependent pathlength changes the solution to the 
propagation equations (|1| and |7|), with the replacements 



where w{T) = (dT/dt)/p IS MV- 

We have allowed the normalization parameter A in A(T) to vary, and 
figures ^ and |8| compare the resulting dLi/ cffie and dB / cffie to the same ratios 
obtained with a constant A(T) = Ao- Note that the changes are relatively 
small. This is because the energy dependence of the pathlength is most im- 
portant at high energies, above the region R(T) ~ A(T) at which the flux 
turns over. However, the majority of the flux resides in the region around 
and below this maximum. Since these energy regimes are largely unaffected, 
the LiBeB production is weakly affected as well. For the momentum source 
spectrum the results with and without an energy-dependent pathlength differ 
by a constant factor of roughly 25%, while for the energy source spectrum the 
difference varies but is always small < 5%. The smallness of the difference 
between the two confinement behaviors justifies the A(T) = A approxima- 
tion used in this work and elsewhere. The systematic decrease of the yields 
for a given scale Aq follows from the energy dependence, under which A has 





(31) 
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a maximum, at T — Ti, of only O.8A0. Thus the confinement is always less 
than this value, and very much lower for the important low energies. 

We have investigated the effect of inelastic nuclear reactions between 
the primary cosmic rays and the ISM on the effective confinement given 
in eq. fl28l) . As discussed by Malaney and Butler (1993), this effect con- 
tributes an extra species-dependent and energy dependent term to the es- 
cape length. Because the pathlength is energy dependent, one must use eq. 
fl3T|) to compute its effect. We have done so, using cross sections given in 
Meyer (1972), Townsend & Wilson (1985), and Bystricky, Lechanoine-Leluc, 
& Lehar (1987). The results appear in Figure [], and the behavior is as ex- 
pected. In the usual case, dLi/dBe decreases with increasing A e ff = A esc , 
and we now have A e ff < A esc . Thus, we expect dLi/dBe to be consistently 
larger at a given A esc , as is the case. This effect is more pronounced at high 
A esc , where the nuclear loss term dominates. 

Another property of the early galaxy that affects confinement is the state 
of ionization of the ISM. Energy losses are increased in an ionized medium, 
leading to a shorter ionization range. Because the confinement always ap- 
pears as R/A (for A(T) constant in energy), this effectively increases the es- 
cape length PCV account for ionization by doubling the energy losses for neu- 
tral media, a procedure which amounts to assuming a fully ionized medium 
(Ginzburg & Syrovatskii 1964). If indeed the ionization state evolves during 
the Pop II epoch, the effective change in A can affect the cosmic ray spectrum 
independent of the considerations in eqns. (|28|) and (|3~T1). 



6.2 Spectrum 

Both WSSOF and PCV have independently argued for a flatter spectrum in 
early cosmic rays than that presently observed. The key observation support- 
ing this argument is the B/Be ratio. Although the observational uncertainty 
here is clearly large, the current preliminary numbers are low compared to 
the present day observed value of 15 ± 3 (with e.g. Mathews, Walker, and 
Viola (1985) predicting ~ 17 for Pop I), and to the Pop II value of ~ 14 
calculated in WS (using a transported spectrum 4>(T) oc (T + m p )~ 2 - 7 and 
Pop I C:N:0 ratios). 

The argument for a flat early cosmic ray spectrum arises if one takes a low 
central value for the Pop II B/Be ratio seriously. (Remember values below 7 
are impossible in any spallation model, and all hypothetical stellar processes 
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suggested thus far act to increase B/Be.) Be and B arise from the same 
spallation processes, so there are a limited number of available differences 
to exploit in fixing the ratio. Their differences in mass and charge lead to 
different ionization losses and ranges, and consequently different confinement 
for the two elements (at a given energy per nucleon). However, with the 
ionization ranges scaling law of eq. fl5|), it is clear that Rb{T) < R Be (T), and 
thus energy loss and thermalization effects (c.f. eq. ||]) favor B rather than 
Be and so act to increase the B/Be ratio. The other important difference 
between Be and B production rates lies in their spallation cross sections. As 
argued in section |6.1| , the differences make the B/Be ratio sensitive to the 
cosmic ray spectrum, with a lower B/Be for flatter spectra. 

Another argument for a flatter spectrum has been offered by PCV (but 
weakened by Malaney and Butler's (1993) arguments; see Prantzos & Casse 
(1993)). PCV note that the low observed Li/Be ratio (once the primordial Li 
component has been subtracted) is more easily reproduced by a flat spectrum 
because of exponential dropoff of a + a cross sections. 

In addition to the effect of changes in the confinement, the propagated 
spectrum obviously depends on the source spectrum. PCV choose as a source 
spectrum a power law in momentum, q(p) oc p~ 2 , a type of spectrum appear- 
ing in models of cosmic ray shock acceleration. Such models assume acceler- 
ation via scattering off of supernova shock waves, and have been successful 
in reproducing galactic and heliospheric spectra in detail (c.f. the Blanford 
& Eichler (1987) review). Note that the PCV spectral index of 2 is the uni- 
versal limiting value for very large Mach number shocks. We note, however, 
that despite the success of and active interest in this model, there remains 
disagreement upon the the appropriate model for acceleration. It is thus 
worthwhile to examine the sensitivity of the LiBeB yields to the assumed 
spectrum, where there remains some flexibility to choose the spectral index. 

Most important for LiBeB production is the choice of spectral type. While 
shock acceleration models suggest source power laws in magnetic rigidity and 
hence momentum are good candidates, the measured (and hence propagated) 
cosmic ray flux at earth is also consistent with a source spectrum in total 
energy per nucleon. The two are of course difficult to distinguish at the rel- 
ativistic energies at which the cosmic ray measurements are unaffected by 
solar modulation effects. Furthermore, solar effects introduce large uncer- 
tainties at lower energies, rendering unreliable the data that could make this 
distinction. In addition, a power law in total energy has the advantage that it 
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remains finite at low energies, thus providing a mockup of physical processes 
responsible for preventing a divergence in the source at low energies. 

Thus we have tried source spectra of power laws both in momentum and 
in total energy per nucleon, i.e. 

q(p) oc p~ 7 

q(T) oc (T + m p )(2m p T + T 2 )-^ +1)/2 , and (32) 
q(T) oc (T + m p )-"' (33) 

with m p the proton mass. For the shock acceleration model-inspired mo- 
mentum spectrum, the index 7 is given by 

■< x { 2 + Tp) (34) 

for the case of large Mach number M = v s hock/c soun d, the ratio of shock ve- 
locity to the sound speed in the shocked gas. Note that the Mach number 
depends upon the local ISM sound speed, which in turn depends upon, e.g. 
the ISM density and pressure. If the early galaxy, or its primordial proto- 
galactic building blocks, were to be significantly hotter or denser than today, 
then potentially the cosmic ray spectrum could be quite different. 

For the total energy spectrum we take 7 = 2.7 , which is consistent 
with the measured, transported cosmic ray spectrum. Note however, that 
this choice is inappropriate for the simple leaky box with energy dependent 
escape length. In this case, we live in the confinement region (the galaxy), 
and so the escape energy dependence is folded into that of the source. The 
proper source index is thus 7 source = lobs — ft ~ 2.1, with (3 the high-energy 
index of the energy-dependent A from eq. (j30l). This is not the case in, 
for example, the nested leaky box model, in which we do not live in the 
confinement region (which shrouds the sources), and so we do not see its 
energy dependence: 7 sowce = jobs- Again, the model uncertainty argues for 
examination of a range of source indices. 

Note the dramatic effect that the spectral index has on the Li/Be and 
B/Be plots (figs. [3] and |j). The large low energy flux of the momentum spec- 
trum makes this spectrum softer than one in total energy, and the resulting 
elemental ratios behave as one would expect for such a soft spectrum. As 
discussed above in section f|, low energy resonance behavior is emphasized 
and leads to large Li/Be and B/Be ratios. 
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We have allowed the spectral index 7 to vary over the range 2 < 7 < 3. 
Results appear in figures 10 and 11 from which it is apparent that the LiBeB 
production rates from a source with a momentum spectrum are more sensitive 
to changes in spectral index than are the rates for a source total energy law. 
Indeed it is clear from fig. ^ that for source spectra in momentum, indices 
with 7 > 2.5 lead to Li-to-Be relative production rates in ratios approaching 
10 3 . This is clearly undesirable given that all observed Be abundances are 
> 1CT 13 , thus implying associated Li production larger than the observed 
plateau, even without the large primordial Li component. 

While the relative B and Be rates (figure [11]) do not vary as dramatically 
for the momentum spectrum, here too they are much more sensitive than the 
same ratios given a total energy spectrum. Additionally, the same range of 
indices that overproduce Li relative to Be give the largest values for B relative 
to Be as compared to low values for the current preliminary but uncertain 
data. Nevertheless, it is intriguing that the range of indices that fail for the 
Li to Be ratio seem to also do the worst for B to Be. Thus momentum spectra 
with indices steeper than about 2.5 are less favored at the present time. Note 
that we get no such constraints for source spectra in total energy, as we see 
only mild sensitivity to a source index for Li to Be (for which all values are 
completely acceptable), and completely insensitivity for B to Be. 

In figure [12] we plot dLi/dBe vs A, as in figure [3], but with different 
spectral indices. For each spectral type, the production ratio grows with the 
steepness of the spectrum, the growth being more pronounced for the scale- 
free momentum law. Figure 13 is a similar plot for the dB/dBe ratio. Note 
that this ratio is mildly sensitive to the momentum spectral index, which 
can change dB/dBe by 1 to 2 units, i.e. 10-20%. Such a change might be 
detectable as observations improve. The dB/dBe ratio is even less sensitive 
to the index of a source with a total energy spectrum. The change over the 
whole range is of order 5%. 

As we have alluded to earlier, a direct comparison of the data and the 
model calculations given here is extremely difficult but can be suggestive. 
The ratio of of Li/ Be given in table 1 contains the primordial Li component 
as well as the cosmic-ray produced Li. In some cases this correction might be 
quite large (a factor of 10 or more) so that the ratios in the table represent 
extreme upper limits. If we had a specific cosmic-ray nucleosynthesis model, 
one could, as shown in Olive and Schramm (1992) extract a Big Bang Li 
abundance. For the cosmic-ray nucleosynthesis model used in WSSOF this 
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led to a value [Li]_B_e = 2.01 ± 0.07, however as we can see from figures 
|5|a and |5|&, there is significant model dependence in the Li/Be ratio. Also, 
the absolute [Li] values are plagued by the poorly given systematic errors 
mentioned in section [|. Furthermore, as we have also noted, the element ratio 
will also be affected by chemical evolution. Thus one can with certainty only 
require the data points as given in the table to lie above any model curve 
in the figures. The trend in the data is also seen in the figure [5], i.e. at 
high metallicities the Li/Be ratio is smaller. Note that a large primordial 
subtraction would favor the lower ratios found with the cosmic-ray energy 
source spectrum. 

To get a feeling for the comparison between data and these models, we 
can make a primordial subtraction corresponding to [Li]^_B = 2.0. It hard to 
make the primordial 7 Li abundance much smaller (Walker et al. 1993; Krauss 
& Romanelli, 1990) and if it were much larger there would be no room for 
cosmic-ray produced Li. In Figure |n] we show the values of Li/ Be with the 
above primordial subtraction. Error bars are too large to display. We also 
plot the ratio from the energy and momentum-spectra for A = 100 gem -2 as 
a function of [Fe]. 

6.3 Other Model Features 

The conventions for reporting spallation/fusion cross sectional data (e.g. Read 
& Viola 1984) are strictly speaking incomplete for calculations of LiBeB pro- 
duction, and thus require a model to make some assumption regarding the 
kinematics of these interactions. Specifically, the tabulated data and semi- 
empirical fits give the total cross section for LiBeB production at a given 
incident energy. As seen explicitly in eq. (Q), however, we require the full, 
differential cross section da^/dT LiBeB (T LiBeB ,T C R). Meneguzzi et al. (1971) 
and many subsequent authors address this issue by assuming that the differ- 
ential cross section is very sharply peaked around those values of the LiBeB 
energy for which T LiBeB = T CR . They put 

da A - 

— {TuBeBiTcRj = a^TcR) 5 (T LiBeB — T C r) ■ (35) 



■ LiBeB 



with afj(TcR) the tabulated cross section. 

While the approximation of sharp peaking over a small range of product 
energies appears to be borne out experimentally in the cases where it has 
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been checked, the other approximation, that of equal energies (per nucleon) 
in the initial and final state is only valid for the case of heavy cosmic rays 
on light targets. We have relaxed the latter assumption, allowing for the 
CR energy to be related to the LiBeB energy via some function Ef(TcR), 
changing eq. ( j35|) to 

da A - 

— ^- (T UBeB , T CR ) = <4(T CR ) S (T LlBeB - E f (T CR )) . (36) 

U-L LiBeB 

With this more general formalism we have investigated two possible kine- 
matic behaviors. One ansatz is that in the center of momentum frame the 
LiBeB is produced at rest, with the lighter debris moving in such a way as 
to conserve energy. This is the assumption we have made throughout this 
paper. One should compare this to another common assumption that in the 
frame in which the light target is at rest, the kinetic energy per nucleon of 
the daughter nucleus is the same as that of the heavy CR parent. We have 
plotted the dLi/dBe ratio for these two cases in figure [15|. 

The issue of spallation kinematics is raised if one considers an alternative 
to the usual production scenario we have presented so far. Duncan et al. 
(1992) and PCV both attempt to address the question of the Be vs Fe slope 
with the suggestion that perhaps most of the LiBeB production is not from 
light (p, a) cosmic rays on heavy (CNO) ISM nuclei, but the reverse: heavy 
cosmic rays on light ISM nuclei. Furthermore, if the heavy cosmic rays do 
not have the ISM abundance or something near to it, then this could yield 
significantly different Be and B production. Increased cosmic ray heavies 
would result, for example, if supernovae were to accelerate their own ejecta. 
The evidence is not definitive on this issue, but there are possible problems 
in maintaining consistency with an assumed first ionization potential depen- 
dence of the cosmic ray source abundances. 

At any rate, we have investigated the dependence of the spallation yields 
on forward and reverse kinematics. Figure [IB] compares the dLi/dBe ratio for 
the cases of normal kinematics (i.e. LiBeB synthesis primarily from light CR 
on heavy ISM with a small admixture of the reverse) to reverse kinematics, 
in which we assume He no ~ Vcnoi an< ^ thus dominate the LiBeB production. 
Note that here the dLi/dBe ratio is much smaller, as the a + a process is not 
as dominant as in the usual early galactic case. 

Note that, in considering the thermalization of LiBeB, authors from 
Meneguzzi et al. (1971) to PCV apparently use S = 1 over all energies 
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in the case of light cosmic rays on heavy targets. These authors note that 
the daughter nuclei move more slowly than the cosmic ray parent and so 
are likely to be stopped. However, as we are considering cases where the 
flux can be very hard, then we may still have quite energetic daughter nuclei 
and it is important to allow properly for incomplete stopping. Note that Li, 
which is made primarily by a + a fusion, is thus the only spallogenic element 
we consider for which (in the dominant production process) the cosmic ray 
projectile not much lighter than the ISM target. Consequently, the Li will 
be the fastest daughter nuclei, and thus accurate treatment of trapping is 
particularly important. To compare with previous work, we have computed 
yields with S = 1 for both kinematics. Results appear in figure [17|. 

We see in figure [T7| that the direction of the effect depends upon the 
spectral type. For a momentum source, one gets a higher dLi/dBe ratio 
when allowing for realistic trapping of LiBeB. The change, ranging from 
about 40% at large A down to 10% at small ones, arises because the Li-to-Be 
production ratio is the smallest in the energy regimes; however reactions at 
these energies produce fast products which are less likely to be thermalized 
than slow ones. Put differently, if Th% ^> Ti ow are different LiBeB energies, 
then in the language of eq. ||, S(Thi) <C S(Ti ow ) This is not the whole story 
for a source law in energy. Here, the spectral peak is much more pronounced 
and at a high energy (~ m p ) than the (lower energy) peak of the flux from 
the momentum source. The flux from the energy source thus rises steeply at 
the (~ lOOMeV) dropoff in the Li production cross sections. Consequently, 
the majority of Li production occurs in the beginning of the high energy 
tail of the cross section. By disfavoring the high energy products through 
incomplete stopping, one removes the faster Li made at the highest overlap 
of cross section and flux. With S — 1, however, these fast daughters are 
included and so raise the Li production relative to Be. 

Let us now look at the dependence of the dLi/dBe production rate on 
the asymptotic value of cross section for a + a — > 6,7 Li. This is a fusion pro- 
cess and is physically different from the other spallation processes. The data 
show that these cross sections decrease exponentially with energy, although 
the highest energy at which a definite cross section has been reported is only 
about 50 MeV/nucleon. Evidence for the continued dropoff comes instead 
from upper bounds to the cross sections, which go out to ~ 200 MeV/nucleon, 
with no plateau detected. In their cross section tabulation, Read and Viola 
(1984) point out the uncertainty in the high-energy a + a reaction, and rec- 
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ommend adopting an exponential decrease with energy, though they do allow 
for the possibility that a plateau might nevertheless eventually be found to 
exist. We feel that an exponential dropoff is the more appropriate descrip- 
tion of a fusion reaction, but some work (notably PCV) has employed such 
a plateau. We have therefore examined the dLi/dBe ratio with and without 
such a plateau, and find that it can have a significant effect, particularly for 
cases of large confinement or spectral index near 2, as these emphasize the 
high energy events in which the differences are important. Figure [18] com- 
pares the dLi/dBe ratio with and without the cross section plateau. While 
the differences are not large for a momentum source (< 10%), for a source 
power law in total energy, the effect is about a factor of 2. Since we see no ev- 
idence as yet for a plateau, in all our other computations none is assumed to 
exist. However, as the effect of such a plateau is significant, we feel it is very 
important for experiments to investigate the behavior of the a + a — > 6 ' 7 Li 
at high energy. 

Finally, we note the important but not decisive contribution to the over- 
all uncertainty of experimental uncertainties in cross sections and in energy 
losses. Cross section uncertainties are typically at a level of 5a /a ~ 10 — 20%. 
While these errors are large enough to be significant, they clearly are smaller 
than the other model effects we have noted above, consequently making 
the model choices important ones. Energy losses are better known, with 
5b jb ~ 5R/R ~ 5% at low energies (where they are unimportant), and < 2% 
at higher energies. These errors do not significantly contribute to our overall 
uncertainty. 

7 Gamma Ray Production 

The same passage of cosmic rays through the ISM that produces LiBeB must 
also create 7-rays, as noted by Silk and Schramm (1993), Fields, Schramm 
and Truran (1993), and Prantzos and Casse (1993). The primary 7-ray source 
from the nuclear cosmic ray component^ arises from inelastic scattering of 
cosmic rays on the ISM, with p + p —>■ tt° — > 27. Thus the 7-ray spectral 

4 As opposed to the cosmic ray electrons, which produce lower energy 7-rays by 
brehmsstrahlung. We do not explicitly consider the electron component here, but Prant- 
zos and Casse argue that it should scale with the nuclear component as it does today, and 
using this scaling they calculate the contribution to the isotropic 7-ray background. 
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peak is around 

The total (isotropic) 7-ray production rate by the nuclear component of 
the cosmic rays, per second per cubic centimeter, is given by 



in an obvious notation. Note that the calculation is similar to that for LiBeB 
production. Indeed this process occurs under very similar conditions to Li 
production, as both processes derive from reactions between protons and a 
nuclei, which in the early galaxy are undepleted both in the cosmic rays and 
in the ISM. Thus we may expect to scale the 7-ray production to that of Li 
via a {a + a — > Li) to the a(p + p — > 7) cross section ratio. One should note, 
however, that such a simple scaling ignores the different energy dependences 
of these processes, and ignores the issue of the incomplete trapping of the Li 
fusion products. 

Because 7-ray production so closely follows that of Li, the model depen- 
dences will be similar. In particular, the effects of cosmic ray evolution will 
be similar, as will the uncertainties in scaling the 7-ray production to the 
Be. In a subsequent paper on specific galactic evolution models we intend to 
quantify the 7-fluxes and use them as a complementary probe to the LiBeB 
abundance behaviors described here. 

8 Conclusions 

In this paper we showed in detail that there are many uncertainties about 
early cosmic ray behavior and therefore a range of possible assumptions from 
which one must choose in building a model for this behavior. Furthermore, 
the variance in possible early cosmic ray models can have a significant effect 
of on the calculation of cosmic ray production of LiBeB (and 7-rays) in the 
early galaxy. The variance in the models derives both from uncertainties 
about present cosmic rays and from uncertainties about the evolution of 
their injection and acceleration, as well as confinement. 

Because these uncertainties are difficult to resolve, we conclude that fu- 
ture work in modeling early cosmic rays must allow for cosmic ray model 
variances and evolution. As shown in section |6.2| , the choice of source spec- 
tral type is significant yet not settled even for today's cosmic rays, with power 




(37) 
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laws both in total energy per nucleon and in momentum being allowed. Also, 
we show that cosmic ray evolution derives not only from changes in source 
strength, but also from changes in confinement, as detailed in sections |] and 



Q| . Evolution of confinement can have important effects, for example in 
altering the naively expected Be-to- Fe slope, as outlined here and shown in 
the detailed model of PCV. Additional considerations are noted in section 



6.3, most notably the effects of inelastic nuclear collisions on confinement, as 



have been noted by Malaney and Butler (1993). 

In light of our findings we now turn to the work to date on early cosmic 
ray production of LiBeB, namely the relatively model-independent work of 
WSSOF and its sequels, and the detailed cosmic ray and chemical evolution 
model of PCV. The flatness of the curves in figure [T3| justifies the claim in 
WSSOF that the B/Be ratio in cosmic ray models is indeed insensitive to the 
details of cosmic ray models. In addition, because Be and B are pure spalla- 
tion (as opposed to a + a fusion) products, their production ratio dB/dBe is 
independent of [Fe/H] and varies only with the C:N:0 ratios of targets. To 
the extent that this ratio remains fairly constant, we can conclude that the 
B/Be ratio is insensitive to chemical evolution as well. We consequently find 
B/Be to be well-chosen as a model-independent signature of the spallation 
process. 

The Li isotopic ratio, discussed in SFOSW, 6 Li/ 7 Li shows a similar inde- 
pendence to cosmic ray model features but does have a moderate sensitivity 
to chemical evolution via the ISM ratio He/CNO, which sets the relative 
contributions of the a + a fusion versus spallation contributions to 6 Li/ 7 Li. 
We find that <i 7 Li/<i 6 Li is always within ~ 30% of unity for all spectra and es- 
cape parameters considered here. This uncertainty is much smaller than that 
of the 6 Li/ 7 Li ratio of Smith et al. (1992). Consequently we find SFOSWs 
assumption of a temporally constant <i 6 Li/<i 7 Li production ratio to be appro- 
priate. 

However, we find that the "zeroth order model" of WS, WSSOF, and 
SFOSW to be overly simplistic not only in its lack of chemical evolution, but 
also in its lack of any kind of cosmic ray evolution. The effect on predictions 
of this model is unimportant for the B/Be limits, but their calculations of the 
Li/Be ratio are, as we have noted, only a lower bound to the actual ratios. 

The PCV model does allow for chemical as well as cosmic ray evolution, 
and in particular they allow for evolution of cosmic ray confinement. We 
thus find their model indeed to be a useful examination of LiBeB production. 
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Our caution is that their model is a specific one made with particular, albeit 
reasonable, assumptions. Their specific cosmic ray source model we find to 
be well-motivated and reasonable, but not a unique choice and we would 
urge that future work consider different momentum source indices, as well as 
source spectra in total energy. Furthermore, PCV admit that their particular 
implementation is to some degree arbitrary. In the first place, it is not 
obvious which parameters control A. Secondly, even given a particular scaling 
prescription, it is not trivial to model the input parameters (e.g. galactic 
scale height and magnetic field) accurately. PCV assume A depends upon 
the scale height of the collapsing disk, but they do not compute the collapse 
explicitly, and so impose an ad hoc time dependence. Clearly a firm model 
for A evolution is crucial, but at present is unavailable. We agree with PCV 
that more theoretical work in this area is needed. 

In summary, we feel that improved data on PopII LiBeB abundances 
with consistent controls on systematics might lead to an ability to constrain 
models for the origin and propagation of cosmic rays in the early galaxy. 
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FIGURE CAPTIONS 



Figure 1: (a) The cosmic ray proton flux and a spectra compared to 
source spectra (q(p) oc p~ 2 ) scaled as Aq. Note that in the limit of no 
energy losses (a good approximation at high energy) we have <fi = Aq, 
and so departure from this scaling indicates the effect of energy losses. 
Note also that the a flux scales very well with the proton flux as the 
He cosmic ray abundance 0.08. That this proportionality is accurate 
over the whole energy range reflects the similar energy loss behavior 
due to the scaling A/Z 2 . 

{b) As in (a), for q(T) oc (T + m p )- 2 - 7 . 

Figure 2: (a) The transported cosmic ray proton flux for a range of confine- 
ment parameters: A = 10, 30, 100, 300, and 1000 gcm~ 2 . The source 
is a power law in momentum, q(p) oc p~ 2 . 

(b) As in (a), for a source power law in total energy, q(T) oc 
(T + m p )~ 2 - 7 . 

Figure 3: The Li to Be production ratio as a function of pathlength. The 
solid line is for a source spectrum that is a power law in momentum, 
q{p) P~ 2 ; the dashed line is for a source power law in total energy, 
q(T) oc (T + m p )- 2 - 7 . We put [C/H] = [N/H] = -2.5, [O/H] = -2.0. 

Figure 4: The B to Be ratio as a function of pathlength, plotted as in fig. 
I 

Figure 5: (a) The production ratio of Li to Be as a function of metallicity 
[Fe/H], for a source q(p) oc p~ 2 , and for a range of conement parameters: 
A = 10, 30, 100, 300, and 1000 gem" 2 . See discussion in the text 
regarding comparison to observations. 

(b) As in (a), for a source q{T) oc (T + m p )~ 2 - 7 . 

Figure 6: A plot of the pathlength distribution of equation fl30|) , chosen in 
Garcia- Munoz et al. (1987) to provide the best fit to the observed B/C 
ratio. 

Figure 7: (a) The dLi/dBe ratio for A(T) = A (solid line), and for A(T) 
given by equation (|30|), with A as indicated on the figure (broken line). 
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The source spectrum is q(p) oc p 2 , and ISM abundances are [C/H] = 
[N/H] = -2.5 and [0/H] = -2.0. 

(b) As in (a), for a source spectrum g(T) oc (T + m p )~ 2 - 7 . 

Figure 8: (a) The cffi/cffie ratio for A(T) = A (solid line), and for A(T) 
given by equation (|30|) , with A as indicated on the figure (broken line). 
The source spectrum is oc p~ 2 and ISM abundances are [C/H] = 
[N/H] = -2.5 and [O/H] = -2.0. 

(b) As in (a), for a source spectrum g(T) oc (T + m p )~ 2 - 7 . 

Figure 9: (a) The dLi/dBe ratio for A = A esc , the usual case (solid line), 
and for A(T) given by eq. (P5|), (broken line), where we have plotted 
against A esc . The source spectrum is q(p) oc p~ 2 and ISM abundances 
are [C/H] = [N/H] = -2.5 and [O/H] = -2.0. 

(6) As in (a), for a source spectrum g(T) oc (T + m p )~ 2 - 7 . 

Figure 10: The ratios of Li to Be production rates, shown as a function 
of source spectral index 7 for a source spectrum in momentum and in 
total energy. Note that the yields from the momentum spectrum are 
very sensitive to the index adopted, and the ratio of rates can become 
unacceptably high in this case. In contrast, the yields from the total 
energy spectrum are relatively insensitive to spectral index, and are all 
at acceptably low levels. 

Figure 11: As in fig. for the ratio of B to Be rates; the same trends 
are seen. Note however that the B/Be ratio is negligibly sensitive to 
the total energy index, and is only mildly sensitive to that for the 
momentum law. 

Figure 12: (a) Curves of dLi/dBe, as a function of A, for different spectral 
indices 7. The source is q{p) oc p~~ 2 . 



(b) As in (a), for q(T) oc (T + m 



1-2.7 



v) 



Figure 13: As in fig. [L2|, for the dB/dBe ratio. 

Figure 14: Similar to figure |5|. We plot dLi/dBe for each source type at 
A = 100gcm~ 2 . The data, as described in the text, has been corrected 
assuming the minimal allowed primordial 7 Li production, [ 7 Li]^_B = 
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2.0 and thus (linearly) subtracting this amount from the Li vales of 
table 1. See discussion in the text regarding caveats in comparing the 
theoretical curves, which are upper bounds, to the data. 

Figure 15: The dLi/dBe ratio for different determinations of the daughter 
nucleus kinematics. The solid curve is calculated assuming that the 
daughter nucleus is stationary in the center of momentum frame of 
the parents; this is the standard used throughout this paper. The 
dashed curve is calculated assuming that, in the rest frame of the light 
parent, the daughter kinetic energy per nucleon is equal to that of the 
heavy parent. The source spectrum is q{p) oc p~ 2 for curves (a and (b, 
and q{T) oc (T + m p )~ 2 - 7 for curves (c and (d. Additionally, we put 
[C/H] = [N/H]=-2.5, [O/H]=-2.0. 

Figure 16: (a) The dLi/dBe ratio for (solid line) cosmic ray CNO abun- 
dances equal to ISM abundances {Ucno = vIjno^ the USU& 1 case), 
and for (dashed line) cosmic ray CNO abundances fixed to be solar 
(Ucno ^ UcMo)- The source spectrum is q(p) oc p~ 2 and ISM abun- 
dances are [C/H] = [N/H] = -2.5, and [O/H] = -2.0 

(b) As in (a), for a source spectrum q(T) oc (T + m p )~ 2 - 7 . 

Figure 17: The dLi/dBe ratio with incomplete LiBeB trapping, i.e. with 
^LiBeB — exp(i?/A) (solid line), and with complete LiBeB trapping 
(•SLiBcB = 1). The source spectra are as labeled, and we put and 
[C/H] = [N/H]=-2.5, [O/H]=-2.0. 

Figure 18: (a) The dLi/dBe ratio in the absence (solid line) and presence 
(broken line) of a high energy plateau in the a + a — > 6 ' 7 Li cross section. 
The source spectrum is q(p) oc p~ 2 and ISM abundances are [C/H] = 
[N/H] = -2.5 and [O/H] = -2.0. Note the small (~ 10%) increase with 
the plateau. 

(b) As in (a), for a source spectrum q(T) oc (T + m p ) -2 ' 7 . Here 
the difference is larger (a factor of ~ 2), as the flux is harder than that 
from the momentum flux, and so the high energy cross section behavior 
is more important. 
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